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摘 要 : 针对 蝴蝶 优化 (monarch butterfly optimization，MBO) 算法 易 陷 入 局 部 最 优 和 收敛 速度 慢 等 问题 ， 提 出 了 
一 种 基于 改进 的 交叉 迁移 和 共享 调整 的 蝴蝶 优化 (MBO with cross migration and sharing adjustment，CSMBO) 算法 。 
首先 ， 利 用 基于 维度 的 垂直 交叉 操作 来 替换 标准 MBO 算法 的 迁移 算 子 ， 形 成 交叉 迁移 算 子 ， 有 效 提升 其 搜索 能 
其 次 ， 将 原始 调整 算 子 改 为 具有 信息 分 享 功 能 的 共享 调整 算 子 ， 以 加 快 算 法 的 收敛 速度 ; 最 后 ， 采 用 贪 焚 选择 策略 
取代 标准 MBO 算法 中 的 精英 保留 策略 ， 减 少 一 次 排序 操作 进而 提高 其 计算 效率 。 为 了 验证 CSMBO 算法 的 优化 能 
力 ， 测 试 了 其 在 30 维和 50 维 的 函数 优化 ， 并 与 三 种 优化 算法 进行 比较 ， 其 实验 结果 表明 CSMBO 算法 具有 良好 的 
优化 性 能 。 
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人 Abstract: In order to solve the problems of the monarch butterfly optimization (MBO) algorithm, such as it is easy to fall 
into local optimum and the convergence Speed is low, this study proposed an improved MBO algorithm with cross migration 
and Sharing adjustment (CSMBO) . Firstly, this paper introduced a dimension-based vertical crossover operation to 
substitute the original migration operation of MBO, and then generated a cross migration operator. Thus, this operation 
could improve search ability of MBO algorithm effectively. Secondly, in order to Speed up the convergence of MBO 
algorithm, the sharing adjustment operator with information sharing replaced the original adjustment operator. Finally, this 
paper utilized the greedy strategy to instead of the elite strategy of MBO, which could reduce one sorting operation and 
improve the calculation efficiency of MBO algorithm. To evaluate the optimization ability of our CSMBO algorithm, this 
paper made some experiments on a set of common benchmark functions with 30-dimensions and 50-dimensions, and the 
results showed that the proposed CSMBO algorithm had good optimization performance, and outperformed currently 
available three optimization approaches with which it is compared. 
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0 引言 optimization, PSO) [2250、 烟 花 算法 031、 生 物 地 理学 优化 算法 
I0、 蝴 蝶 优 化 〈(MBO ) 算法 05、 莉 火 虫 优化 算法 0618 等 。 

优化 是 在 某 些 约束 条 件 下 搜索 和 找到 给 定 问题 的 最 佳 解 而 在 这 其 中 ， 蝴 蝶 优 化 算法 因 其 优化 性 能 优异 而 得 到 广泛 关 

的 过 程 趾 。 对 于 优化 算法 ， 它 们 可 以 分 为 确定 性 算法 和 随机 注 。 
算法 两 大 类 中。 近年 来 ， 许 多 可 以 转换 成 优化 问题 的 现实 问 蝴蝶 优化 算法 是 2015 年 由 Wang 等 人 05 受 自然 界 黑 脉 金 
题 变 得 越 来 越 复杂 , 并 且 很 难 通 过 确定 性 算法 来 解决 B-53。 因 降 蝶 迁徙 行为 启发 而 提出 的 。MBO 算法 有 两 个 重要 的 算 子 
比 ， 在 实际 应 用 中 对 随机 算法 的 研究 已 成 为 一 个 很 有 前 景 的 即 迁移 算 子 和 调整 算 子 。 前 者 提供 一 定 的 局 部 搜索 能 力 ， 后 
热点 课题 1。 如今, 群 智能 优化 算法 作为 随机 算法 的 一 个 重 者 提供 一 定 的 全 局 搜索 能 力 。 在 MBO 算法 中 ， 蝴 蝶 的 搜索 
要 分 支 , 已 展示 出 解决 复杂 优化 问题 的 强大 而 有 效 的 能 力 氏 方向 由 迁移 算 子 和 调整 算 子 确定 ， 而 这 两 个 算 子 可 以 同时 执 
其 中 比较 著名 的 有 : 粒子 群 优化 算法 (particle swarm 行 ， 因 此 ，MBO 算法 很 适合 用 于 并 行 处 理 ， 并 且 能 够 很 好 
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松 症 的 分 类 问题 请 ! 等 


仍 存在 搜索 全 2 较 权 和 


于 MBO 人 数 少 且 易 于 
很 多 领域 ， 如 0-1 背包 问题 09、 自 
、 动 态 路 径 问 题 久 、 骨 质 下 G 
但 MBO 算法 在 解决 一 些 复杂 问题 时 


因此 ， 许 多 学 者 对 其 进行 了 改 


进 研 究 。Wang 等 人 [9 提 


基于 自 适应 交叉 算 子 和 贪 
禁 策 略 的 MBO (monarch butterfly optimization with greedy 


strategy and self-adaptive crossover operator, GCMBO) 算法 ， 


易 陷 入 局 部 最 优 和 收敛 速度 
林 优 化 问题 时 ， 以 各 奉 仿 大 的 下 迁 守 
本 文 针 对 MBO 算法 易 陷入 局 
题 ,提出 了 一 种 基于 改进 的 交 
(CSMBO) 算法 。 首 先 ， 


提高 了 MBO 算法 的 优化 能 力 ; Feng 等 
斯 变异 引入 MBO 算法 ， 增 强 了 MBO 算法 的 全 局 优化 性 能 ; 
蒙 丽 萍 等 人 29 采用 将 种 群 动态 随机 分 
各 ， 以 保持 种 群 搜索 的 多 样 性 ; 
的 MBO 算法 ， 以 帮助 MBO 算法 跳 
等 人 2 将 反 向 学 习 和 高 斯 扰动 引入 
MBO 算法 的 收敛 速度 ， Hu 等 人 [3 
各 的 MBO 算法 ， 提 高 了 MBO 算 
MBO 算法 也 经 常 和 WA 


是 出 了 一 种 


人 的 将 混沌 理论 和 高 


将 MBO 算法 和 人 工蜂 群 优 


吉 合 。Ghanem 等 人 B0 


法 结合 ， 以 解决 数值 优化 问 
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MBO 算法 的 两 个 位 置 更 新 策略 如 下 : 


割 成 两 个 子 群 体 的 策 
Feng 等 人 提出 了 一 种 多 策略 
出 局 部 最 优点 ?1; Feng 
MBO 算法 ， 从 而 加 快 
自 适 应 种 群 策 
法 的 探索 能 力 等 。 此 外 ， 


:2 ee 在 1+1 时 刻 的 第 K 维 ; 
i 在 上 时 刻 的 第 上 维 。 当 前 迭代 次 数 
区 红 al 和 r2 分 别 从 Subpopulation1 和 Subpopulation2 中 随机 
Ph peri 为 迁移 期 ， 在 MBO 算法 中 取 
peri= 1.2，rand 是 [0, 1] 中 的 随机 数 。 


。 r=randxperi ， 其 吕 


1 NP2 = NP— NP!1, 


7 
| 如 erS 忆 
ME 

Xp 人 > 已 


题 ，Ghetas 等 人 5 将 和 声 搜索 算法 与 MBO 算法 进行 结合 ， 
是 高 了 MBO 算法 的 优化 性 能 。 


上 述 一 些 算 法 仍 存在 


共享 


用 基于 维度 的 垂直 交叉 操作 蔡 换 原始 迁移 操作 ， 形 成 新 的 交 
又 迁移 算 子 ， 以 增加 种 群 多 样 1 
以 在 一 定 程度 上 避免 MBO 算法 陷 


将 原始 调整 算 子 改 为 具 
效 地 利用 种 群 的 有 效 信息 ， 


其 是 在 解决 一 些 复 


部 最 优 和 收敛 速度 慢 等 问 
调整 的 新 的 MBO 
针对 MBO 算法 中 的 迁移 算 子 ， 运 


洋 账 * 


生 ， 提 升 其 搜索 能 力 ， 这 样 可 
局 部 最 优 的 问题 ; 其 次 ， 


ee 
Xjk = Xpest,k 


条 件 下 ,， 若 rand >BAR， 则 


了 信 ， 息 分 享 功能 的 共享 调整 算 子 ， 有 
ti MRO 算法 的 食 敛 速度 ; 


最 后 ， 采 用 贪 禁 选择 策 
策略 ， 减 少 一 次 排序 操 
率 。 为 了 验证 CSMBO 久 法 的 


住 PWo 算法 中 的 精英 保留 


基准 函数 上 进行 测试 实验 ， 
复杂 问题 上 具有 良好 的 优化 性 能 
1 标准 MBO 算法 

MBO 的 相关 知识 可 以 参见 文献 [15，24]。 
中 ， 所 有 的 蝴蝶 个 体 都 是 到 


大 南部 《区 域 1)、 墨 西 哥 〈 


为 调整 率 ， 在 MBO 算法 
于 随机 数 rand, 在 1+1 时刻 的 第 天 维 就 会 更 新 ,其 


高 了 MBO 算法 的 计算 效 
i 化 能 力 ， 在 一 组 不 同 维度 的 
其 结果 表明 CSMBO 算法 在 解决 


其 次 ， 其 他 蝴蝶 的 位 置 通 


过 调整 算 子 来 更 新 。 所 以 ， 央 


蝶 个 体 的 搜索 方向 是 通过 迁移 


算 子 和 调整 算 子 确定 的 ， 而 且 迁 
执行 。 因 此 ，MBO 算法 非常 适合 并 行 处 理 ， 能 很 好 地 
局 探索 能 力 。MBO 算法 遵循 如 下 理 


想 规则 : 


平衡 局 部 搜索 能 力 和 全 


a) 所 有 的 蝴蝶 都 只 
蝴蝶 组 成 了 整个 蝴蝶 种 群 。 


b) 每 一 只 蝴蝶 后 代 都 是 区 志 


子 产生 的 。 


c) 为 了 保持 种 群 数 量 不 变 ， 
有 相应 的 父 代 蝴蝶 消失 。 

d) 具有 最 好 的 适应 度 值 的 蝴蝶 
对 其 进行 任何 操作 。 这 样 可 以 保证 
迭代 次 数 的 增加 而 下 降 。 


2 中 的 蝴蝶 通过 迁移 入 


且 有 后 代 蝴 蝶 产 生 ， 


体 


子 和 调整 算 子 能 够 同时 


子 a 的 计算 方式 为 a= 呈 


在 MBO 算法 
E 想 化 的 ， 仪 位 于 美国 北部 和 加 拿 
区 域 2) 两 个 区 域 。 因 此 ， 蝴 蝶 
的 位 置 以 迁移 算 子 和 调整 算 子 两 种 方式 更 新 。 首 先 ， 后 代 的 
产生 (位 置 更 新 ) 通过 迁移 算 子 ; 


从 式 (3) 可 以 看 出 , 第 + + 1 


x+ 一 0.5) 
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的 是 更 新 蝴蝶 在 区 域 1 和 2 之 间 的 迁 
先 ， 蝴 蝶 总 数 为 VP， 区 域 1 和 2 
NP1 = ceil(p x NP) 和 和 


迁移 率 )，MBO 算法 设置 p = 方 


或 等 于 x 的 最 接近 的 整数 。 
区 域 2 的 子 群 记 为 Subpopulation2。 因此 ， 


FP 的 蝴蝶 数 分 别 为 
其 中 jp 为 贿 蝶 比 


ceil() 是 舍 入 到 大 于 


区 域 1 的 子 群 记 为 Subpopulation1 ， 
迁移 算 子 可 以 表 


(1) 


V4 表示 碌 在 1 时 刻 的 
1 表示。 


的 是 更 新 Subpopulation2 中 蝴蝶 的 位 


二 于 蝴蝶 j 的 所 有 维 , 若 rand <p， 则 该 个 体 进行 位 置 更 


(2) 


PP; 对 为 加 在 1+ 1 时刻 的 第 维 ; Wwx 为 Landl 与 Land2 


; er 为 Xbest 的 第 天 维 ; 1 洲 当 前 迭代 次 

则 该 个 体 进行 位 置 更 新 的 公式 如 下 : 
X= (3) 
为 的 第 天 维 ， 是 从 Land2 中 随 出 。 在 此 


， 于 是 可 知 这 种 


原 迁 移 算 子 有 机 


尺 蝴蝶 的 信息 
直接 继承 的 方法 过 于 简单 ， 没 
华 的 有 效 信息 。 因 此 ， 将 共享 操作 与 原 调整 算 子 结 


融合 


中 权重 
，Smax 为 蝴蝶 的 最 大 步 长 ，dx 为 
过 Levy flight 计算 dx = Levy( *)。 

基于 交叉 迁移 和 共享 调整 的 蝴蝶 优化 算法 
改进 算法 的 动机 

MBO 算法 的 迁移 算 子 采 用 


日 形 


， 提 升 算 


新 的 共享 调整 算 子 ， 这 样 可 以 有 效 地 利 月 


在 区 域 1 和 2 中 。 即 区 域 1 和 2 中 的 


该 个 体 进行 位 置 更 新 的 公式 如 下 : 


(4) 


Ph BAR = p。 如 果 BAR 


呀 


式 (1) 所 示 的 直接 取代 式 迁移 
0 定 的 局 部 搜索 能 力 ， 但 迁移 方式 单一 ， 
局限 ， 故 将 交叉 操作 与 
成 交叉 寺 萝 算 子 ， 这 样 可 以 有 效 增加 种 群 的 多 样 1 


,是 直接 从 第 4 
有 充 


种 群 的 优 


于 引导 蝴 


全 MBO 算法 的 精英 保留 策 
曾 加 了 算法 计算 的 复杂 程度 。 因 此 ， 月 


ny 
HT 


就 会 


自动 进入 下 一 代 , 不 
j 蝶 种 群 的 质量 不 会 随 着 


度 的 垂直 交叉 操作 取代 


各 需要 设置 参数 也 


以 减少 一 次 排序 操作 ， 进 而 降低 算法 的 计算 复杂 度 。 
2.2 ”交叉 迁移 算 子 

鉴于 MBO 算法 的 迁移 算 子 搜索 能 力 不 足 ， 运 用 
其 迁移 算 子 的 直接 取代 式 迁 移 操作 ， 
形成 新 的 交叉 迁移 算 子 ， 可 以 有 效 弥 补 原始 迁移 操作 迁移 下 


换 精英 保留 策略 ， 不 仅 可 以 避免 精英 参数 的 设置 ， 而 


蝶 向 最 优 解 的 方向 靠近 ， 进 而 加 快 算 


上 排序 操 
有 贪 禁 选择 策略 蔡 
有 还 可 


基于 维 


-WW 


一 以 及 在 解 空间 中 可 搜索 到 的 位 置 较为 


局 
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的 不 足 。 下 面 给 出 新 的 交叉 迁移 算 子 ， 其 计算 公式 可 以 表示 
为 : 

加 Cxxipt+(l—-Oxxiygrsp 

2 ge G5) 


其 中 : 垂直 交叉 缩放 因子 c=rand, q=randxD。 
对 式 (1) 和 式 (5) 的 分 析 可 知 ， 在 式 (5) 中 ， 第 t+ 1 代 蝴 蝶 
的 第 维 会 受到 第 t 代 蝴蝶 的 第 Ek 维和 第 gq 维 (9 是 随机 选择 
的 ) 两 者 的 共同 影响 ， 其 具体 方式 是 对 两 者 进行 加 权 交 叉 计 
算 , 这 样 可 以 增加 种 群 的 多 样 性 , 进而 增强 算法 的 搜索 能 力 ， 
有 效 克 服 式 (1) 存 在 的 形式 简单 、 方 向 单一 等 缺陷 。 
另外 , 由 于 贪 禁 选择 策略 具有 只 接受 更 优秀 个 体 的 优点 ， 
接 下 来 可 以 将 原 迁 移 算 子 中 的 精英 保留 策略 更 换 为 贫 禁 选择 
策略 ， 进 而 降低 交叉 迁移 算 子 的 计算 复杂 度 。 其 表示 公式 如 
下 : 


(6) 


其 中 : xiww 是 下 一 代 新 生成 的 蝴蝶 个 体 ，fGm) 和 f(D 分 别 
表示 蝴蝶 x 和 xi 的 适应 度 值 。 
F 面 给 出 交叉 迁移 算 子 的 伪 代 码 ， 如 算法 1 所 示 。 
算法 1 交叉 迁移 算 子 
fori=1toNPl do 
fork=1toDdo 
根据 式 (5) 更 新 砍 
end for k 
根据 式 (6) 计 算 xi 
end fori 
2.3 共享 调整 算 子 
针对 调整 算 子 中 式 (3) 直 接 继承 方式 过 于 简单 的 缺点 ， 为 


new 


了 ps 
x = 


Xt!,Otherwise 
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式 如 下 : 

Xi = Xhosx + (rand —0.5) x (Xhosx 一 2 
其 中 :xo 为 种 群 中 最 好 的 个 体 ws 在 上 时 刻 的 
是 [0, 1] 中 的 随机 数 。 

从 式 (7) 可 以 看 出 ,共享 调整 算 子 充分 利用 种 群 的 有 效 信 
息 ， 有 利于 引导 蝴蝶 向 最 优 解 的 方向 靠近 ， 加 快 算法 的 收敛 
速度 ， 同 时 ， 也 将 原始 调整 算 子 中 的 精英 保留 策略 改 为 贪 焚 
选择 方案 来 提高 共享 调整 算 子 的 计算 效率 ， 其 表示 公式 如 式 
(6) 所 示 。 

下 面 给 出 共享 调整 算 子 的 伪 代 码 ， 如 算法 2 所 示 。 
算法 2 共享 调整 算 子 
forj= lto NP2 do 
for 上 = ltoDdo 
取 一 个 [0, 1] 区 间 的 随机 数 rand 
if rand <p then 
根据 式 (2) 计 算 zx 
else 
从 子 群 
( 记 作 3) 
根据 式 (3) 计 算 鸡 
if rand >BAR then 
根据 式 (4) 计 算 区 
end 让 
end 让 
end for k 
根据 式 (6) 计 算 ev 
end forj 
2.4 CSMBO 算法 总 流程 


(7) 


rand 


江 


大 维 ; 


2 (Subpopulation2) 中 随机 选 出 一 只 蝴蝶 


了 充分 利用 种 群 中 的 有 效 信 息 , 可 以 将 共享 操作 合并 到 MBO 在 提出 的 交叉 迁移 算 子 和 共享 调整 算 子 的 基础 上 ， 给 出 
算法 的 调整 算 子 中 ， 进 而 形成 新 的 共享 调整 算 子 。 其 表示 公 CSMBO 算法 的 流程 如 图 1 所 示 。 
Ca 
设置 允 歼 并 初始 化 种 尖 
令 1=1 
(< 
机 
根据 每 只 蝴蝶 的 位 置 计算 各 自 的 适应 度 什 
vv 
| 根据 蝴蝶 的 适应 度 值 进行 排 ] 
机 
把 蝴蝶 种 群 分 成 两 个 子 群 


人 


根据 算法 1 更 新 Subpopulation1 


根据 算法 2 更 新 Subpopulation2 
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| t=t+1 | 
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<1<MaxGen? > 


图 1 CSMBO 算法 流程 


Fig.1 Main procedure of CSMBO algorithm 
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3 ”仿真 实验 可 以 增加 寻 优 难度 ， 以 便 验 证 CSMBO 算法 处 理 各 种 复杂 优 
化 问题 的 能 力 。 四 种 算法 的 共同 参数 设置 如 表 2 所 示 。 其 中 ， 

3.1 测试 函数 与 参数 设置 最 大 迭代 次 数 MaxGen 随 优化 函数 的 维 数 增加 而 增加 ， 其 他 
为 了 测试 本 文 提出 的 CSMBO 算法 的 有 效 性 ， 在 一 组 不 参数 设置 均 与 相关 文献 的 PSO、MBO、GCMBO 保持 一 致 。 

同 维度 的 基准 函数 上 进行 四 种 算法 的 优化 实验 对 比 。 表 1 中 在 表 2 中 ， 选 择 两 个 维度 (30 维和 50 维 ) 进行 算法 性 能 测 
给 出 了 八 个 测试 函数 的 信息 。 它们 包括 单 峰 函数 如 fi~ 户 , 多  ” 试 ， 其 中 30 维 对 应 的 种 群 数量 和 最 大 友 代 次 数 分 别 是 50 和 
峰 函 数 如 f~ 记 。 其 中 ， 单 峰 函 数 只 有 一 个 全 局 最 优点 ,适合 1 000，50 维 对 应 的 种 群 数量 和 最 大 迭代 次 数 分 别 是 50 和 
用 于 评估 算法 的 搜索 能 力 ; 多 峰 函 数 有 多 个 局 部 最 优点 ， 适 ”2 000。 在 实验 实施 过 程 中 ， 若 种 群 规模 过 小 ， 会 导致 算法 收 
合用 于 评估 算法 的 探索 能 力 。Fmin 是 函数 的 理论 最 优 值 。 本 ”全 速度 过 慢 ;， 若 种 群 规模 过 大 ， 会 增加 算法 的 时 间 复 杂 度 ， 
文 所 有 仿真 实验 均 在 操作 系统 为 Windows 7、CPU 为 主 频 。” 浪费 不 必要 的 存储 空间 。 在 文献 [24] 提 出 的 结合 贪 禁 策略 和 
3.10 Ghz 和 内 存 为 4GB 的 PC 上 实施 ， 编 程 语言 采用 自 适应 交叉 算 子 的 MBO 算法 、 文 献 [26] 提 出 的 采用 动态 分 割 


MATLAB R2014a。 种 群 策 略 的 改进 MBO 算法 以 及 文献 [29] 提 出 的 基于 自 适 应 
在 实验 中 ， 将 本 文 提出 的 CSMBO 算法 与 其 他 三 种 相关 种群 的 改进 MBO 算法 中 ， 均 将 种 群 规模 设置 为 S0， 因 此 ， 
的 优化 算法 (PSO 、MBO03、GCMBOBP9) 进行 优化 性 能 本文 在 分 配种 群 规模 时 ， 考 虑 到 基准 函数 的 维度 和 MBO 算 
对 比分 析 。 本 文 用 四 种 算法 对 不 同 维度 的 基准 函数 进行 优化 ”法 的 特点 ， 也 将 种 群 规 模 设 置 为 S0。 为 了 避免 随机 性 对 实验 
实验 ， 以 便 验 证 CSMBO 算法 的 有 效 性 、 高 优化 效率 和 普 适 结果 的 影响 ， 在 仿真 实验 中 ， 四 种 算法 对 每 一 个 优化 问题 均 
性 等 。 函 数 的 维度 选择 包括 30 维和 50 维 。 不 同 的 维度 选择 ” 独立 进行 30 次 实验 。 
表 1 基准 函数 信息 
Table 1 Information of benchmark functions 
函数 公式 定义 域 in 
Quartic f = Vi +random[0,1) [~—1.28,1.28]? 0 
Schwefel 1.2 fo (7) -sy [=100, 100]? 0 
Schwefel 2.21 户 CO =maxi(x|l,1<i<d) [=100, 100]? 0 
Step fi(W)= Ys +0.5))? [=100,100]? 0 
Rosenbrock fs(#) Ss D? +100G% —xXn)?) [-10, 10]? 0 
fi =F00sini ry) + Sy —D*(1+10sin? (zy)) + (ya 一 D2) 
+ Du,10,100,4), 
Penalized 1 部 [=50, 50]? 0 
k(x —a)”,x >a, 
其 中 ，y -1 了 Gi +1), u(xi,a,k,m)=10,-a< x <a, 
Kk(—x; —a)”, xX <—a, 
(x)=0.1(sin? (Axx ) 十 $0 —1)?(l+sin?(3xx,)) 
Penalized 2 . [=50, 50]? 0 
(x —D(+sin?(2xx,)) + Pu,5,100, 4) 
Griewank 及 CD 2 ep [-600, 600]? 0 
表 2 四 种 算法 的 参数 设置 值 上 是 最 好 的 ， 但 是 CSMBO 算法 的 最 差 值 、 平 均值 和 标准 
Table 2 Values of parameters of four algorithms 差 在 四 种 算法 中 都 是 最 优 的 。 值 得 一 提 的 是 ， 在 方 和 foc 上， 
种 群 数量 维度 最 大 迭代 次 数 CSMBO 算法 的 平均 值 和 标准 差 明显 比 其 他 三 种 算法 的 结果 
50 30 1000 更 好 。 在 及 和 有 上，CSMBO 算法 的 优化 结果 均 比 其 他 三 种 
50 50 2000 算法 优秀 ， 特 别 是 在 户 上 ，CSMBO 算法 的 结果 大 幅 领 先 于 
3.2 优化 性 能 比较 其 他 三 种 算法 。 在 及 和 包 上 ，MBO 算法 、GCMBO 算法 和 
在 8 个 基准 函数 上 对 不 同 维度 (30 维和 50 维 ) 的 测试 ”CSMBO 算法 在 最 优 值 的 寻 优 上 均 取 得 了 一 致 的 结果 ， 尤 其 
结果 如 表 3 和 4 所 示 。 其 中 包含 30 次 独立 运行 寻 优 结果 中 的 是 在 户 上 ， 都 找到 了 函数 的 理论 最 优 值 ， 但 是 CSMBO 算法 
最 优 值 、 最 差 值 、 平 均值 和 标准 差 ， 粗 体 表示 测试 结果 的 最 的 最 差 值 、 平 均值 和 标准 差 的 结果 都 是 最 好 的 。 总 之 ， 从 30 
优 者 。 文 献 [9]、 文 献 [15] 和 文献 [24] 均 指出 平均 值 和 标准 差 维 函 数 上 的 实验 结果 可 知 ， 本 文 提 出 的 CSMBO 算法 的 优化 
值 越 小 说 明 算 法 的 搜索 能 力 越 好 且 稳 定性 越 强 。 性 能 相对 优异 。 
表 3 给 出 了 四 种 算法 在 30 维 函数 上 的 测试 结果 。 从 表 3 表 4 给 出 了 四 种 算法 在 50 维 函数 上 的 测试 结果 。 从 表 4 
中 的 数值 可 以 看 出 ， 在 户 和 上方 ~ 户 上 ，GCMBO 算法 在 最 优 中 可 以 看 出 ， 对 于 有 和，CSMBO 算法 的 寻 优 结果 均 是 最 


优 的 。 在 户 上 ，PSO 算法 的 最 差 值 和 标准 差 是 最 好 的 ， 而 
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CSMBO 算法 的 最 优 值 和 平均 值 是 最 优 的 。MBO 算法 、 总 之 ， 从 整体 情况 来 看 ，CSMBO 算法 在 八 个 测试 函数 上 的 
GCMBO 算法 和 CSMBO 算法 均 找 到 了 户 的 理论 最 优 值 ， 但 收敛 性 能 都 是 最 优 的 ， 再 一 次 证 明 CSMBO 算法 的 改进 策略 
CSMBO 算法 的 最 差 值 、 平 均值 和 标准 差 结 果 都 大 幅 优 于 其 ” 是 有 效 的 。 其 原因 是 由 于 CSMBO 算法 融入 了 交叉 操作 ， 增 
他 三 种 算法 。 对 于 广 、 训 和 方 , GCMBO 算法 的 最 优 值 结果 是 。 强 了 种 群 的 多 样 性 ， 提 升 了 算法 的 全 局 勘探 能 力 ; 而 且 在 调 
最 好 的 ， 而 CSMBO 算法 在 最 差 值 、 平 均值 和 标准 差 的 评价 。” 整 算 子 中 加 入 了 共享 操作 ， 种 群 中 更 优 个 体 的 信息 得 到 有 效 


指标 上 是 最 优秀 的 。 值 得 一 提 的 是 , 在 户 和 所 上 ,CSMBO ”利用 ,使 得 CSMBO 算法 的 收敛 速度 得 到 了 很 大 幅度 的 提升 。 
算法 的 结果 都 大 幅度 优 于 其 他 三 种 算法 。 总 之 , 从 50 维 函 数 综 上 所 述 ，CSMBO 算法 不 管 在 单 模 态 函 数 上 ， 还 是 在 
上 的 实验 结果 可 知 ， 本 文 提出 的 CSMBO 算法 的 优化 性 能 是 。 多 模 态 函数 上 ,与 其 他 三 种 算法 相 比 具有 更 优秀 的 优化 性 能 ， 
相对 优秀 的 。 所 以 本 文 CSMBO 算法 的 改进 是 有 效 和 可 行 的 。 
表 3 四 种 算法 在 30 维 函 数 上 的 测试 结果 表 4 四 种 算法 在 50 维 函 数 上 的 测试 结果 
Table 3 Results of four algorithms on 30-dimensional functions Table 4 Results of four algorithms on 50-dimensional functions 
函数 ”算法 最 优 值 最 差 值 平均 值 标准 差 函数 ”算法 最 优 值 最 差 值 平均 值 标准 凌 
PSO 1.6029 12.6050 6.9126 2.5498 PSO 21.1397 57.5429 39.4603 8.8743 
MBO 0.0015 122.1438 50.5840 45.3491 MBO 0.0034 415.5476 153.3220 142.7727 
7 GCMBO 3.3127e-24 7.7247 1.5967 2.2949 7 GCMBO 1.917le-07 133.2460 19.750 29.8177 
CSMBO 4.5748e-20 1.6384 0.0629 0.2998 CSMBO 1.89Sse-17 0.6991 0.0437 0.1377 
PSO 8.8927e+03 1.8163e+04 1.3409e+04 2.5485et+03 PSO 2.4468e+04 4.7567e+04 3.6691e+04 5.6002e+03 
MBO 7.1783 6.5083e+04 2.4090e+04 1.8531e+04 MBO 0.0010 1.6716e+05 7.1776e+04 5.1240e+04 
7 GCMBO 3.8183 2.6657e+04 1.1791e+04 8.5563et+03 7 GCMBO 2.3034e+03 7.4799e+04 2.8123e+04 1.7575e+04 
CSMBO 4.3512e-08 0.0015 6.3360e-05 2.7170e-04 CSMBO 1.2244e-07 5.7843e+04 S$.2787e+03 1.5786e+04 
PSO 37.4256 54.7068 46.3699 4.5268 PSO 56.5356 94.9000 87.4604 11.0036 
MBO 0.0010 87.6419 33.0509 22.3718 MBO 0.0071 90.5485 35.3305 26.7642 
h GCMBO 0.0947 73.3843 31.8467 23.9338 记 GCMBO 0.0299 85.3000 33.9794 21.5405 
CSMBO 6.6164e-05 12.2039 2.4087 3.5827 CSMBO 5.241le-05 15.5230 4.1900 4.9262 
PSO 10379 16015 1.3538e+04 1.4832e+03 PSO 26696 33563 3.0441e+04 1.9358e+03 
MBO 0 72991 2.5781le+04 2.7767e+04 MBO 0 132209 5.1815e+04 5.1859e+04 
f GCMBO 0 2206 203.7667 541.4205 A GCMBO 0 12346 1696 2.9548e+03 
CSMBO 0 人 0.2333 1.2780 CSMBO 0 19 0.8333 3.6016 
PSO 2.3207e+t06 2.3380e+07 9.0459e+06 5.2705e+06 PSO 1.9350e+07 8.9184e+07 4.9437e+07 1.8547e+07 
MBO 1.2839e-1ll 5.8334e+08 1.3545e+08 1.9159e+08 MBO 3.4589e-ll 1.114le+09 2.4720e+08 4.0804e+08 
GCMBO 2.5213e-16 8.2390e+06 5.6040e+05 2.0417e+06 GCMBO 2.0184e-17 5.4740e+07 3.2179e+06 1.0925e+07 
CSMBO 6.9181e-12 11.6346 0.8624 2.3041 CSMBO 1.0028e-11 41.S839 4.7557 9.4187 
PSO 7.5324e+06 8.1079e+07 4.2706e+07 1.9612e+07 PSO 8.6540e+07 2.4986e+08 1.5685e+08 4.5239e+07 
\ MBO 9.8230e-10 1.0516e+09 1.4534e+08 2.6536e+08 MBO 6.5719e-10 1.9497e+09 3.6l3le+08 6.4416e+08 
/ GCMBO 4.7593e-14 5.2262e+07 3.4594e+06 1.1075e+07 / GCMBO 9.8724e-16 1.4513e+08 1.0977e+07 3.1121le+07 
CSMBO 8.1086e-11 1.2936 0.0914 0.2833 CSMBO 2.3708e-10 194.3978 10.7529 37.6554 
PSO 412.4324 1.5276e+03 963.1479 230.7740 PSO 1.3892e+03 3.5147e+03 2.4754e+03 478.7173 
MBO 3.9119e-04 9.5470e+03 1.8755e+03 3.1023e+03 MBO 1.3455e-04 1.8360e+04 5.9326e+03 6.2480e+03 
7 GCMBO 1.4106e-10 626.1261 98.7457 142.3047 7 GCMBO 3.392Se-09 1.488le+03 243.2459 353.8553 
CSMBO 1.7249e-05 140.8104 52.2129 38.1776 CSMBO 1.4367e-05 265.6663 65.0644 58.8563 
PSO 91.0693 142.6368 114.4050 12.9686 PSO 200.9274 306.9407 266.5297 26.2755 
MBO 1.0000 595.4902 148.0689 176.6548 MBO 1.2936 1.0722e+03 419.0091 381.3554 
n GCMBO 1.0000 130.8612 42.5119 45.7053 , GCMBO 1.0000 312.7389 96.2010 97.1711 
CSMBO 1.0000 1.0000 1.0000 2.5766e-11 CSMBO 1.0000 1.0000 1.0000 2.2857e-11 
以 上 在 不 同 维度 基准 函数 上 的 实验 结果 证 明了 CSMBO ”3.4 CSMBO 算法 的 复杂 度 分 析 
算法 的 有 效 性 与 可 行 性 ， 不 仅 在 不 同 维度 单 峰 函数 上 有 较 好 在 软 硬 件 运行 环境 相同 的 前 提 下 ， 优 化 算法 的 复杂 度 计 
的 收敛 精度 ， 而 且 在 不 同 维度 多 峰 函 数 上 显示 了 更 好 的 全 局 。 算 主要 由 目标 函数 的 复杂 度 和 算法 自身 流程 的 复杂 度 两 部 分 
搜索 能 力 ， 以 相对 明显 的 优势 优 于 其 他 三 种 对 比 算法 。 组 成 。 本 文 在 对 比 实验 中 ， 设 置 所 有 算法 (PSO、MBO、 
3.3 收敛 性 分 析 GCMBO 和 本 文 的 CSMBO 算法 ) 具有 相同 的 种 群 数量 和 最 


为 了 更 直观 地 看 出 CSMBO 算法 的 收敛 速度 、 局 部 搜索 ”大 迭代 次 数 ， 使 得 它们 的 最 大 函数 评价 次 数 C(MNEFE) B2 近 
能 力 和 全 局 搜索 能 力 , 本 文 给 出 了 四 种 算法 在 求解 50 维 测试 ”” 似 相等 ， 故 CSMBO 算法 的 复杂 度 计算 主要 取决 于 自身 流程 
函数 时 的 收敛 曲线 图 ， 如 图 2 所 示 。 的 复杂 程度 。 假 设 CSMBO 算法 的 最 大 迭代 次 数 为 7， 种 群 
从 图 2 可 以 看 出 ，CSMBO 算法 在 对 八 个 测试 函数 进行 。 规模 为 N, 子 群 1 为 YI， 子 群 2 为 N - Nl1， 维 度 为 D。 由 
优化 时 的 收敛 速度 均 优 于 其 他 三 种 对 比 算法 。 特 别 是 在 广 、 CSMBO 算法 流程 (图 1) 可知, 该 算法 的 时 间 复 杂 度 主要 1 
记 \ 广 和 上 , CSMBO 算法 的 收敛 速度 大 幅度 优 于 对 比 算法 。 ”每 次 迭代 循环 决定 ， 详 细 步 又 分 析 如 下 : a) 需 要 计算 蝴蝶 的 
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适应 度 值 ,时 间 复 杂 度 为 O(N);b) 排 序 , 时 间 复 杂 

时 间 复 杂 度 是 O(N); dd) 运 
时 间 复 杂 度 为 O(N1 x 站 )， 
个 内 循环 ， 时 间 复 杂 


c) 分 成 两 个 子 群 ， 
有 两 个 内 循环 ， 
子 中 ， 也 存在 


度 为 O((N -NN1) x DD)。 


林 ， 等 : 


度 是 OOV2); 
运行 交叉 迁移 算 子 ， 
量 在 共享 调整 算 


间 受 到 种 群 


所 以 , CSMBO 算法 的 总 时 间 复 杂 度 是 TCD = OUUD) = OCT x 
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(N+N+N+(NIxD)+((N-NI)xD)))= OTxQN+N+ 
NxD)) = O(Tx NN”)。 在 CSMBO 算法 中 ,由 
规模 N 和 变量 维度 D 的 影响 , 则 
度 为 S(n) = O(ftn)) = O(N x D)。 


于 可 变 的 储存 空 
计算 其 空间 复杂 
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Fig. 2 Convergence curves of four algorithms on eight benchmark functions 
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本 文 针 对 MBO 算法 易 陷入 


题 ， 提 


局 部 最 优 和 收敛 速度 慢 等 问 
出 了 一 种 基于 交叉 迁移 和 共 
(CSMBO) 算法 。 首 先 ， 将 MBO 的 迁移 算 子 改 为 基于 旨 


< 时 
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